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ABSTRACT 

In the course of the Sloan Digital Sky Survey (SDSS-I), a large fraction of the surveyed 
area was observed more than once due to field tiling overlap, usually at different epochs. We 
utilize some of these data to perform a supernova (SN) survey at a mean redshift of z = 0.2. 
Our archival search, in ~ 5% of the SDSS-I overlap area, produces 29 SN candidates clearly 
associated with host galaxies. Using the Bayesian photometric classification algorithm of Poz- 
nanski et al., and correcting for classification bias, we find 17 of the 29 candidates are likely 
Type la SNe. Accounting for the detection efficiency of the survey and for host extinction, this 
implies a Type la SN rate of r Ia = (14.0+^+} '\ ± 2.5) x 10~ 14 h% yr" 1 L~] g , where the errors 
are Poisson error, systematic detection efficiency error, and systematic classification error, re- 
spectively. The volumetric rate is R ia = (l.S^o^aii ± 0A2 ) x 10 ~ 5 y f_1 h io M PC~ 3 . Our 
measurement is consistent with other rate measurements at low redshift. An order of magni- 
tude increase in the number of SNe is possible by analyzing the full SDSS-I database. 

Key words: Supernovae: general Cosmology: observations, miscellaneous Surveys 



1 INTRODUCTION 

Supernovae (SNe) play a central role in galaxy evolution and cos- 
mic metal production. Measuring the rates at which SNe explode is 
thus an important step for understanding the chemical evolution of 
the universe. In recent years, efforts have intensified to measure the 
low-redshift Type la SN rate both in field environments (Cappellaro 
et al. 1999; Hardin et al. 2000; Madgwick et al. 2003; Blanc et al. 
2004, Botticella et al. 2008, Dilday et al. 2008) and in galaxy clus- 
ters (Gal- Yam et al. 2002; Maoz & Gal- Yam 2004; Sharon et al. 
2007; Mannucci et al. 2008; Graham et al. 2008; Sand et al. 2008). 
However, due to small SN numbers, there are still significant un- 
certainties in low-redshift SN rates. 

In this paper, we demonstrate that a large archival repository of 
SNe, one that is potentially useful for a low-redshift rate measure- 
ment using a large number of SNe, exists in the data from the first 
phase of the Sloan Digital Sky Survey (SDSS; York et al. 2000). We 
use a small fraction of these data to detect and compile a sample of 
SNe la and to derive the SN la rate at low-redshift. The techniques 
we use in this paper are also of relevance for future projects such 
as the Panoramic Survey Telescope and Rapid Response System 
(Pan-STARRS; Kaiser 2004) and the Large Synoptic Survey Tele- 
scope (LSST; Tyson 2002). These projects will survey huge areas 
in a relatively short time, and will produce large samples of SNe for 



which spectral classification will not be possible, due to their large 
numbers. 

In §2 we describe the SDSS data we use. The pipeline used to 
process these data and detect SNe is presented in §3. Detection ef- 
ficiency and photometric calibration are discussed in sections 4 and 
5. §6 presents our preliminary results, including a first SN sample, 
its classification, and a calculation of the SN la rate. We compare 
our results to previous measurements in §7, and summarize in §8. 



2 SDSS IMAGING DATA 

The SDSS imaged about one quarter of the sky in five bands (u, g, r, 
U z, centered at 355lA, 4686 A, 6165 A, 7481 A, 8931 A; Fukugita 
et al. 1996). Images were photometrically (Tucker et al. 2006) and 
astrometrically (Pier 2003) calibrated by the SDSS pipeline (Lup- 
ton et al. 2001). The data products of the SDSS (images and object 
catalogs) were made availabl^ in a series of Data Releases (see 
Adelman-McCarthy et al. 2008 for a description of the latest data 
release, DR6). 

With the objective of covering the survey area once, imaging 
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Figure 1. Examples of overlap between SDSS-I fields centered around 
RA= 14:56:49, Dec=+1 0:57:48. The dashed polygon is an overlap between 
two adjacent strips which is used for SDSS quality checks and the solid 
polygon is an overlap resulting from the mapping of the celestial sphere on 
to rectangles. 

was performed by scanning the sky in great circles. Each scan was 
along a 2. 5 -degree- wide strip, where each strip was divided into nu- 
merous "fields". However, dividing the celestial sphere on to rect- 
angular planes causes the rectangles to overlap, especially close to 
the poles of the survey scan coordinate system (see e.g. Fig. 1). In 
addition, adjacent strips have an intentional overlap for the purpose 
of photometric and astrometric quality checks. The fact that differ- 
ent strips were imaged at different times raises the possibility of 
using the overlap regions to detect transient events. 

Each final SDSS field is an image of 2048 x 1361 pixels. The 
image pixel scale is 0"396 with a median point spread function 
(PSF) of 1"4 in the r band. An exposure time of 53.9 sec was 
used to image all fields, resulting in point source AB magnitude 
95% repeatability limits of 22.0, 22.2, 22.2, 21.3 and 20.5, in the 
w, g, r, /, z bands, respectively. 



3 SUPERNOVA SURVEY DATA PIPELINE 

To deal with the vast amounts of data in the SDSS database, we 
developed a largely automatic pipeline for downloading individual 
subsets of overlapping field images, and processing them one at a 
time. Our pipeline consists of three independent modules for down- 
load, registration, and detection, executed in that order. 

As a first step, we compiled a list of overlapping SDSS fields. 
To do so, we downloaded the coordinates of all the fields in the 
SDSS DR4 database. By applying a polygon intersection algo- 
rithm, which assumes planar geometry, to the list of coordinates, 
we constructed a list of the overlapping regions of each SDSS field. 
Each image set, consisting of a first-epoch image (the "reference" 
epoch) and its overlapping second epoch images, was individually 
downloaded for further processing by our pipeline. In the present 
paper, we search for SNe in the region 220° < RA < 240°, and 



— 1° < Dec < 64°. This region is not far from the pole of the SDSS 
coordinate system (RA = 275°, Dec = 0°), resulting in a large over- 
lap area of 92 deg 2 , obtained from 460 deg 2 of SDSS images. 

Our survey search method is based on image subtraction. We 
note that an alternative method is to search for SNe in the SDSS 
catalog using different criteria, e.g., SN colors (see Poznanski et 
al. 2002). However, a SDSS catalog SN search has some disadvan- 
tages. For example, a blind color search will be affected by color 
contamination originating from SN host galaxies. In addition, lack- 
ing a direct access to the SDSS pipeline makes it difficult to esti- 
mate the survey detection efficiency function. 

We chose to limit our SN search to the g and r bands since 
they are the deepest bands in the SDSS. Furthermore, the scanning 
order of each field in the SDSS is r, /, u, z, g. Therefore, the r- and 
g-band exposures of the same field in a given scan have the largest 
time separation, i.e., there is a ~ 5 minute difference between expo- 
sures of the same field in these bands. This time difference is critical 
for identifying and excluding solar- system objects from among the 
SN candidates. 

Two computers were used for running our pipeline. One com- 
puter was used for continuous downloading of images from the 
SDSS database. In parallel, the registration and detection modules 
(see below) were run on a computer with a Pentium IV 3.4 GHz 
processor and 2 GB of memory. The download rate and the pro- 
cessing rate both dictated a net data flow rate of about 1 deg 2 per 
day. In practice, software, hardware, and communications problems 
resulted in a lower rate, and guided our decision to stop the current 
search after about 90 deg 2 . 

3.1 Image registration 

The registration module aligns the overlapping images in each set 
to their reference image and produces a difference image in which 
SN candidates are searched for by the detection module. Both the 
g and r image sets, once downloaded, are registered separately by 
the registration module. An overlapping image is first registered to 
its reference image, based on their world coordinate system (WCS) 
coordinates, using the wregister IRAF (Tody 1986) task. The over- 
lapping areas are then cut out from both the reference and registered 
images. Using each field's photometric parameters, extracted from 
the SDSS databasfl we next linearly match the background level 
and the zeropoint of the registered image to those of the reference 
image. At this point, we try to achieve a better image registration, 
which is crucial for image subtraction, by matching the positions 
of objects which appear in both images. We first detect objects in 
the overlapping image segments by applying the Source Extractor 
(SExtractor) program (Bertin & Arnout 1996) to both segments. 
Next, by cross correlating the object positions, a more precise reg- 
istration between the two segments is obtained using the geomap 
and geotran IRAF tasks, allowing for offsets in the x and y axes 
and a rotation angle between the two images. In order to avoid poor 
statistics in the matching process, this latter alignment is performed 
only if there are at least seven matching objects. Otherwise this 
stage is skipped. 

Next, the image with the smaller PSF full width at half max- 
imum (FWHM) is degraded by convolving it with a 2D Gaussian 
kernel, G(x,y) exp[— (x 2 + _y 2 )/2(7 2 ], in order to match the PSF 
of the second image. The kernel is found from the parameters of 
the two image PSFs, listed as "psfWidth" in the SDSS catalogue 

2 http://casjobs.sdss.org/casjobs 
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"Field" Table. This simplistic PSF matching approach is dictated 
by the small number of objects (generally not point sources) in the 
overlapping regions, which prevents the application of more so- 
phisticated PSF-matching algorithms (e.g. Alard & Lupton 1998; 
Alard 2000). 

Following subtraction of two registered images, the absolute 
values of the difference image is formed, so that all residuals are 
positive. In order to smooth out residuals due to imperfect align- 
ment, the difference image is smoothed by convolving it with a 2D 
Gaussian, three pixels wide (la). 

3.2 SN candidate detection 

The residuals are detected in the difference image by applying SEx- 
tractor to the image. Since the final difference image is positive def- 
inite (see §3.1), it has a one sided noise distribution. We chose to 
apply a 6a detection threshold in the detection process. The value 
of the detection threshold is calculated using the Poisson fluctua- 
tions of the background counts in both the reference and registered 
images. The residuals detected by SExtractor are automatically ex- 
amined in more detail to screen for various non-SN detection con- 
taminations, as described next. 

We first search for variable stars within our candidate list. Us- 
ing SExtractor, we obtain a list of objects in both the registered 
and reference images. If an object is detected in both images at the 
same position where a residual was detected in the difference im- 
age, we query the SDSS catalogue for objects at that position. If 
an object, catalogued as a star, exists at that position, the candi- 
date is considered a variable star and is rejected from our candidate 
list. Similarly, we reject candidates spectroscopically identified as 
quasars. 

We next explore the possibility that a residual is the result of 
poor image registration. We search, using SExtractor, for positive 
residuals in two new difference images: the reference minus regis- 
tered image, and vice versa. If in each of the two images a residual 
is detected near the position of a candidate, we compare the differ- 
ence between the photon counts of the two residuals to our detec- 
tion threshold. In contrast to our original detection in the absolute 
value of the difference image, we now require the photon count dif- 
ference to have at least 3a significance. In a final test for improper 
alignment, a stamp of 41 x 41 pixels around the position of each 
residual is cut out of the registered and reference images. The two 
image stamps are re-registered using the xregister task in IRAF, 
based on cross correlation. A new difference image is produced us- 
ing the new re-registered images. If no residuals are detected in the 
new subtracted stamp image, the candidate is discarded. 

The remaining candidates are subjected to another test, aimed 
at determining whether or not they are moving objects. We first 
query the SDSS catalogue to check if the target has been flagged as 
a moving object. We also compare the position of each candidate 
in the g and r bands, assuming that it was detected in both bands. If 
the candidate position has changed by more than 2 pixels, it is also 
considered a moving object and is excluded from the candidate list. 
All stages up to this point are performed automatically, with no 
human intervention. 

The remaining candidates are saved, together with their sub- 
traction images, for visual inspection, performed by a single per- 
son (AH). The inspection helps reject false positives of various 
types, such as artifacts and residuals due to poor PSF matching, 
poor image alignment, cosmic rays, and saturated objects (see Fig. 
2). About 99% of the candidates found by the automatic pipeline 
are discarded as false positives in the visual inspection stage. 




Figure 2. Examples of false positives which were rejected in the visual 
inspection stage. The left panel shows a residual in the subtraction image 
due to image misalignment. The residual in the right panel is due to poor 
PSF matching. 



4 DETECTION EFFICIENCY FUNCTIONS 

Estimating the efficiency of our SN detection process is critical for 
deriving a reliable SN rate. We have planted a sample of fake SNe 
in the SDSS images, whose recovery fraction provides an estimate 
of the detection efficiency as a function of SN magnitude. 

The first step in producing the fake SN sample was choosing 
the SN hosts. We compiled the g magnitude and the photometric 
redshift (Csabai et al. 2003; Oyaizu et al. 2008) of each galaxy in 
the overlapping sets of images from the SDSS catalog. Under a 
simplifying assumption that the SN rate is proportional to stellar g- 
band luminosity, we selected a random subset of galaxies weighted 
by luminosity. To each of these selected hosts we then assigned a 
fake SN. The SN was assigned a random g-band absolute magni- 
tude in the range of —19.5 to —7.5. The absolute magnitudes were 
converted to observed magnitudes using a distance modulus based 
on the SN host photometric redshift (assuming a Hubble parameter 
of H = 70kms" 1 Mpc- 1 , a mass density in units of the critical 
density £l m = 0.3, and cosmological constant O.^ = 0-7). The SN 
r-band magnitude was randomly chosen to be in the range of —0.5 
and +1.5 of the g band magnitude, a range motivated by calculat- 
ing synthetic g and r magnitudes from a set of observed spectra of 
SNe la (Nugent 2002; Poznanski et al. 2002; Poznanski, Maoz, & 
Gal- Yam 2007). 

The fake SNe were added to the real images as part of the 
data processing, prior to image registration, as follows. First, we 
randomly chose the image in which the SN was to be planted, i.e., 
either the reference image or the registered one. We then cut out 
a region around the fake SN host of size 1.5 times the host's 90% 
light radius. SExtractor was applied to the host stamp image, pro- 
ducing a list of 10%, 20%, 100% light radii of the host. The 
radial distance of the fake SN from its host center was chosen ran- 
domly from among these annuli, assuring that the locations of the 
artificial SNe roughly follow their galaxy host light. The final po- 
sition of the fake SN with respect to the host was at a randomly 
chosen position angle. The SN was then planted in the selected im- 
age using the IRAF task mkob jects. 

The fake SN sample underwent the same processing as the 
real data, including the visual inspection stage, ensuring it reflected 
faithfully the actual detection efficiency. By spreading the fake SNe 
among all the overlapping fields, we also took into account the fact 
that the efficiency may vary from field to field. 

Our detection efficiency functions in the g and r bands are 
shown in Figure 3. We find that our efficiency level is ~ 60%, at 
best. This is probably due to the poor quality of image subtraction 
when the PSF matching and/or the registration are not perfect. For 
example, bright hosts often leave large residuals at their centers in 
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Figure 3. Detection efficiency functions in the g (top) and r (bottom) bands. 
Error bars represent Id Poisson Errors. 



the difference images. Therefore, a real SN that is close to its host 
center, may be mistaken for an artificial residual due to poor reg- 
istration, by either the pipeline or the human inspector. This effect 
seems to be independent of magnitude, even for bright SNe, due to 
the fact that such SNe will tend to be hosted by nearby, and hence 
bright, galaxies. It is also evident that the detection efficiency in 
the r band starts declining at brighter magnitudes, compared to the 
g band, probably due to the relative faintness of galaxies in the g- 
band. 



5 PHOTOMETRIC CALIBRATION 

The SN candidates which pass the visual inspection are further ex- 
plored. In order to obtain the magnitudes of these candidates, we 
download the images, in which a candidate resides, in the remain- 
ing (w, /, z) SDSS bands. Since the images of a field in different 
bands do not fully overlap, we first align them according to the 
image in which the candidate was detected (in either the g or the 
r band). We perform this alignment for the reference field images 
and for the registered field images separately. Then, the reference 
and registered field images in each band are processed by the same 
registration module used to originally process the g and r bands 
(see §3). After a final difference image is obtained in all five bands, 
41 x 41 -pixel images are cut out around the candidate. 
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These images are used to perform aperture photometry of 
the candidates. The counts are summed in an aperture of radius 
2a, where a is 1/2.35 of the FWHM of the average PSF of a 
field (taken from the SDSS catalog "Field" table). With the zero- 
point, airmass, and extinction parameters, also listed in the SDSS 
catalog, we convert the candidate counts to magnitudes. To each 
magnitude we then apply an aperture correction, that accounts for 
the flux outside the aperture. From aperture photometry on bright 
and isolated SDSS stars, we find mean aperture corrections of 
0.21,0.22,0.23,0.25,0.25 in the w, g, r, i, z bands, respectively. 
With these corrections, our final magnitudes for bright stars also 
match those in the SDSS catalogue. We corrected the candidate 
magnitudes for Galactic extinction according to Schlegel et al. 
(1998). 

To obtain realistic error estimates for our magnitude measure- 
ments, we have performed the same analysis on a large sample of 
artificial SNe. The artificial SNe were blindly planted in several 
images in the same manner as done for obtaining our detection ef- 
ficiency functions (see §4). We expect the photometric errors to be 
larger than the usual Poisson errors due to inaccurate registration, 
varying backgrounds, and residuals from the host galaxy subtrac- 
tion. We planted ~ 400 artificial SNe in each of eight magnitude 
bins (see Table 1). In each bin, and for each band, we calculate 
the root-mean square (rms) of the difference between the measured 
magnitudes and the original magnitudes assigned to the artificial 
SNe, and adopt it as the systematic photometric error. The results 
are listed in Table 1 . 



6 RESULTS 
6.1 SN sample 

Our final candidate list for the 92 deg 2 of overlap area searched in- 
cludes 47 transient candidates which we denote as SISN01 to 47, 
where SISN stands for SDSS-I SN. Among the candidates, 25 are 
clearly associated with a detected host galaxy, but offset from the 
nucleus (if there is one). A further 11 are, to within SDSS resolu- 
tion, at the centers of their hosts. A final 1 1 candidates are "host- 
less", i.e., cannot be unambiguously associated with any detected 
galaxy. Our criterion for hostlessness is being separated by both 
> 5" and more than two times the 90% light radius from any galaxy. 

We first turn our attention to estimate the sample contamina- 
tion by SN "impostors". Based on the SDSS limiting magnitudes 
(see section §2), we estimate that a true SN la, at maximum light, 
will be undetectable at redshifts z > 0.35. Although there are ex- 
ceptions, most core-collapse SNe are less luminous than SNe-Ia. 
Candidates with spectroscopic or photometric host redshifts with a 
lcr lower limit above z > 0.35 are therefore excluded. Six of the 
seven candidates excluded by this criteria are at their host centers, 
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and are thus likely to be active galactic nuclei (AGN), rather than 
SNe. The seventh candidate, SISN47, which is not at the center of 
its host, might be a valid SN candidate falsely rejected due to an 
error in the photometric redshift of its host. However, if it were 
included in our final sample it would have no effect on the total 
number of Type la SNe, being classified as a likely core-collapse 
SN (see Table 2). 

The hostless candidates could be either real SNe which reside 
in galaxies below the SDSS limiting magnitude, or they can be im- 
postors such as quasars, slow-moving asteroids, and variable stars. 
To estimate the expected fraction of SNe hosted by galaxies fainter 
than the SDSS limiting magnitude, we use the Blanton et al. (2003) 
galaxy luminosity function in the r band. The fraction of the stellar 
luminosity in galaxies with a luminosity L < L\ im is 

^lim 

/ L(j)(L)dL 

P(L<L lim ) = 4 , (1) 

fL(l)(L)dL 
o 

where is the luminosity function. At z = 0.2, the mean redshift of 
SNe probed by our search (see §6.3, below), the SDSS r-band flux 
limit correspond to an absolute magnitude of M r = — 17.8 mag, and 
P(L < Ln m ) = 23%. Again assuming that SNe track the stellar lu- 
minosity, we therefore expect 23% x 40 « 9 candidates in galaxies 
below the SDSS limiting magnitude, consistent with the 1 1 hostless 
candidates we find. Conversely, this also argues that most of the 
hostless candidates are likely real SNe, as otherwise a large deficit 
of SNe in low-luminosity galaxies would be implied. Nonetheless, 
due to the lack of redshifts for the hostless candidates, we are un- 
able to determine with great confidence which of those candidates 
are real SNe and what are their types. We therefore exclude the 1 1 
hostless candidates from our sample for the purpose of the SN rate 
calculation. The exclusion of the hostless candidates is accounted 
for in the luminosity-normalized SN rate calculation by using the 
luminosity density which originates from galaxies above the SDSS 
limiting magnitude (see §6.3). 

Kuiper belt objects (KBOs) might also play a role as SN im- 
postors. However, they are normally found near the ecliptic, while 
all our SN candidate have ecliptic latitude j5 > 15°. Moreover, ac- 
cording to their magnitude distribution (Bernstein et al. 2004) their 
expected number in our survey, even near the ecliptic, is at least an 
order of magnitude lower than the observed number of SN candi- 
dates. We have checked for asteroids near the positions of all can- 
didates in the Minor Planet Centei0 and Nasa Jet Propulsion Labo- 
ratory databases. No known asteroids were found within 5 r of any 
candidates. Our final main SN sample thus consists of 29 candi- 
dates with hosts, among which five are nuclear and therefore may 
be AGNs. Figures 4, 5, 6 show sections of the reference, registered, 
and difference images for each candidate in the final, hostless and 
high-z samples, respectively. 

6.2 Supernova classification 

Spectroscopic classification of our SN candidates is, of course, im- 
possible in this archival survey. We therefore adopt a photometric 
classification method, the Supernova Automated Bayesian Classifi- 
cation (SN-ABC) routine of Poznanski, Maoz, & Gal- Yam (2007a). 
This method compares the SN candidate magnitudes to a sample of 

3 http://www.cfa.harvard.edu/iau/mpc.html 

4 http://ssd.jpl.nasa.gov/sbfind.cgi 




Figure 4. Final sample of SN candidates - For each candidate the reference 
(left), registered (center), and difference (right) images are shown. Images 
are 16" on a side. Arrows indicate north, with east to the left when facing 
north. 

SN spectral templates of different types, ages, redshifts, and ex- 
tinctions and returns the probability of a candidate being a Type 
la SN, P(Ia), as opposed to being a core-collapse SN. The routine 
uses as a prior the host redshift probability distribution function of 
each candidate. Six candidates in our sample have spectroscopic 
host redshifts. For the rest of the candidates we retrieve photomet- 
ric redshifts from the SDSS catalog. We then assume Gaussian red- 
shift probability distribution functions using the host spectroscopic 
or photometric redshifts and their errors. 

We next perform an analysis similar to that of Poznanski et 
al. (2007b) in order to determine the degree of classification un- 
certainty and bias, i.e., what is the most probable real range in the 
fraction of SNe la, given the fraction that is classified as such. Us- 
ing the SN spectral templates, we create a sample of fake SNe of 
the four most prevalent types: la, IIP, Iln, and Ib/c. The distribution 
of redshifts is according to the redshift distribution of the galaxy 
sample used in the detection efficiency simulations (see §4). We 
note that while Type la's and Iln's are observable up to a redshift 
of z ~ 0.35, given the SDSS limiting magnitude, types Ib/c and IIP 
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Figure 4. [continued] 
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Figure 5. Same as Fig. 4 but for the hostless candidate sample. 




Figure 6. Same as Fig. 4 but for the high-z candidate sample, i.e. the seven 
candidates rejected based on their high photometric host redshifts. 

are observable only out to z ~ 0.15. We find that, in our redshift 
range, the value of P(Ia) is not indicative of the actual probabil- 
ity of an object being a SN la. This is due to the fact that while 
SNe la usually are assigned very high values (85% get P(Ia) > 0.9) 
for core-collapse SNe the results are less obvious. About half the 
type Iln and Ib/c SNe, and a fifth of the type IIP SNe get misclas si- 
fted, with values of P(Ia) > 0.5. As a consequence we use P(Ia) as 
a quality indicator, rather than a probability, cutting the sample at 
P(Ia) > 0.9 where the contamination by core-collapse SNe is min- 
imal. The fractions of simulated SNe that are classified or misclas- 
sified as Type la by the SN-ABC routine using the above threshold 
are 0.85, 0.04, 0.28, and 0. 13 for types la, IIP, Iln, and Ib/c, respec- 
tively. Applying SN-ABC together with a P(Ia) > 0.9 threshold on 
the 29 candidates in our sample results in 16 candidates which are 
classified as la's, with % 2 values in the range 0.3 — 10.6. Table 2 
lists the candidate properties, including SN types and host redshifts. 

Since we do not know the real fraction of Type la SNe at 
low redshift, we create a sample of possible fractions of the dif- 
ferent types, e.g., 30% la, 50% IIP, 10% Iln, 10% Ib/c, using 
steps of 10%. For each set of fractions, we calculate the binomial 
probability of finding the observed 16 SNe la out of 29, using the 
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Table 2. SN candidate sample 
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E{B — V) 


Keasniit 


P(Ia) 


Sample 


SISNOl 


15:44:10.57 


+51:46:30.57 


67.7 


22.43 


20.72 


20.91 


21.04 


— 


0.014 


— 


0.98 


hostless 


SISN02 


15:45:04.35 


+40:11:38.27 


57.7 


20.94 


20.95 


21.54 


21.80 


21.65 


0.016 


0.13 


0.02 


main 


SISN03 


14:56:43.00 


+42:29:17.73 


55.5 


22.40 


21.47 


21.49 


21.22 


21.59 


0.019 


0.20 


0.98 


main 


SISN04 


14:57:06.28 


+46:37:25.51 


59.0 


21.61 


21.31 


21.64 


21.07 


22.38 


0.015 


0.61 


0.42 


high-z 


SISN05 


14:56:48.96 


+10:57:47.53 


26.5 


22.84 


21.18 


21.24 


21.18 


21.41 


0.034 


0.16 


0.99 


main 


SISN06 


14:56:49.05 


+36:36:56.41 


50.3 


22.05 


20.36 


20.56 


20.74 


21.65 


0.015 


0.20 


1.00 


main 


SISN07 


14:42:59.88 


+62:42:48.38 


69.1 


22.44 


20.59 


20.39 


20.69 


20.65 


0.015 


0.16 


0.98 


main 


SISN08 


14:44:51.65 


+06:55:48.55 


21.8 


20.36 


20.39 


20.50 


20.71 


21.25 


0.033 


0.76 


0.01 


high-z 


SISN09 


14:46:19.33 


+53:47:14.39 


63.7 


22.66 


21.04 


20.87 


21.60 


21.22 


0.010 


0.65 


0.93 


high-z 


SISNIO 


14:49:18.82 


+52:33:07.73 


63.0 


23.66 


21.54 


21.35 


21.66 


21.23 


0.015 


0.07 


0.02 




SISN11 


15:00:35.10 


+00:48:53.50 


17.2 


23.63 


20.91 


20.97 


20.74 


21.69 


0.049 


0.20 


0.99 


main 


SISN12 


15:00:58.34 


+29:04:11.96 


43.8 


22.57 


20.98 


20.50 


20.89 


21.25 


0.018 


0.04 


0.19 


main 


SISN13 


15:02:35.07 


+28:56:59.42 


43.9 


21.79 


19.72 


19.30 


19.79 


20.07 


0.023 


0.14 


1.00 


main 


SISN14 


15:03:47.16 


+04:56:04.76 


21.3 


— 


21.96 


20.97 


21.19 


20.82 


0.040 


0.44 


0.39 


main 


SISN15 


15:07:33.29 


+45:33:49.61 


59.1 


22.21 


21.75 


21.31 


21.29 


21.41 


0.021 


— 


0.23 


hostless 


SISN16 


15:08:45.09 


-00:38:53.27 


16.3 


23.41 


21.07 


20.51 


20.33 


21.06 


0.065 


— 


0.74 


hostless 


SISN17 


15:08:37.26 


+45:29:16.72 


59.2 


22.98 


22.20 


21.31 


21.06 


21.31 


0.025 


0.09 


0.93 


main 


SISN18 


15:12:23.54 


+02:40:28.03 


19.8 


22.44 


21.18 


21.15 


21.01 


— 


0.041 


0.23 


1.00 




SISN19 


15:12:44.37 


+06:20:37.23 


23.3 


21.95 


21.52 


20.97 


20.65 


20.17 


0.034 


— 


0.57 


hostless 


SISN20 


15:14:43.66 


+04:39:55.11 


21.9 


— 


21.31 


21.12 


20.91 


20.81 


0.051 


0.10 


0.19 




SISN21 


15:16:38.95 


+05:57:26.38 


23.2 


20.82 


20.13 


19.90 


19.78 


19.87 


0.042 


— 


0.35 


hostless 


SISN22 


15:16:56.21 


+47:10:04.82 


61.4 


21.10 


21.47 


20.86 


21.22 


22.26 


0.030 


0.94 


0.00 


high-z 


SISN23 


15:17:32.35 


+04:30:02.51 


21.9 


21.22 


19.43 


19.31 


19.55 


18.72 


0.047 


0.00 


0.08 




SISN24 


15:17:33.01 


+39:19:52.54 


54.6 


22.24 


20.81 


20.71 


21.09 


21.36 


0.017 


— 


0.95 


hostless 


SISN25 


15:22:29.67 


+38:46:35.4 


54.5 


22.67 


21.68 


21.56 


21.36 


21.72 


0.017 


0.25 


0.99 


main 


SISN26 


15:33:57.56 


-00:48:42.1 


17.9 


21.50 


20.27 


20.25 


20.48 


21.16 


0.116 


0.12 


1.00 


main 


SISN27 


15:33:07.47 


+29:54:54.7 


47.2 


23.18 


22.15 


21.03 


21.22 


21.19 


0.03 


0.08 


0.95 


main 


SISN28 


15:33:40.44 


+09:16:37.1 


27.6 


20.62 


20.57 


20.72 


20.59 


20.96 


0.041 


0.05 


0.02 


main 


SISN29 


15:37:37.48 


-00:38:37.3 


18.2 


22.25 


22.32 


21.02 


20.90 


20.97 


0.098 


0.16 


1.00 


main 


SISN30 


15:43:11.48 


-00:23:55.6 


18.8 


21.70 


20.70 


20.38 


20.14 


19.96 


0.096 


0.04 


0.04 


main 


SISN31 


15:45:27.58 


+26:27:58.8 


44.9 


19.88 


18.61 


18.07 


18.16 


17.82 


0.049 


0.03 


0.01 


main 


SISN32 


15:45:46.02 


+35:37:06.7 


53.4 


21.73 


20.02 


19.55 


19.98 


20.59 


0.029 


0.06 


0.97 


main 


SISN33 


15:46:03.12 


+22:58:53.0 


41.6 


23.01 


22.47 


21.17 


21.21 


20.94 


0.055 


0.12 


1.00 




SISN34 


15:46:15.29 


-00:37:05.8 


18.8 


21.76 


19.83 


19.25 


18.91 


18.94 


0.103 




0.06 


hostless 


SISN35 


15:46:48.36 


+03:26:30.7 


22.8 


20.68 


20.30 


20.83 


20.76 


20.89 


0.094 


0.10 


0.83 


main 


SISN36 


15:49:35.40 


+39:59:10.9 


57.8 


22.28 


21.17 


20.49 


20.15 


20.55 


0.012 


0.56 


0.86 


main 


SISN37 


15:51:22.86 


+04:19:46.6 


23.9 


21.33 


20.71 


20.65 


20.89 


22.35 


0.079 




0.96 


hostless 


SISN38 


15:52:32.30 


+25:38:47.8 


44.6 


22.58 


20.58 


20.02 


20.35 


21.72 


0.060 




0.99 


hostless 


SISN39 


15:52:55.62 


+03:40:15.3 


23.3 


23.30 


20.55 


20.27 


20.62 




0.153 


1.18 


0.97 


high-z 


SISN40 


15:54:12.31 


+24:15:30.8 


43.3 


23.52 


21.85 


21.11 


21.24 


22.42 


0.048 


0.08 


0.59 


main 


SISN41 


15:54:52.01 


+21:07:10.8 


40.4 


24.77 


19.03 


17.90 


17.75 


17.48 


0.055 


0.04 


0.00 


main 


SISN42 


15:54:26.70 


+03:41:09.8 


23.4 


21.85 


20.00 


19.65 


19.57 


19.69 


0.154 




0.64 


hostless 


SISN43 


15:55:53.55 


+31:23:24.4 


50.3 


21.44 


21.08 


21.21 


21.16 


21.10 


0.025 




0.34 


hostless 


SISN44 


15:58:38.83 


+05:15:48.3 


25.2 


22.12 


19.95 


19.22 


19.84 


20.05 


0.056 


0.07 


1.00 


main 


SISN45 


15:59:11.36 


+46:17:49.3 


64.3 


21.66 


19.75 


19.39 


19.97 


20.94 


0.016 


0.15 


1.00 


main 


SISN46 


15:43:30.40 


-01:11:51.3 


18.1 


21.64 


21.00 


21.52 


21.85 


21.05 


0.117 


0.56 


0.79 


high-z 


SISN47 


15:39:04.47 


+03:48:51.4 


22.7 


21.73 


21.22 


20.98 


21.04 


20.20 


0.061 


0.51 


0.04 


high-z 



Notes: Candidates with names in boldface are classified as Type la in the final sample. Magnitudes are before correction for Galactic extinction, derived from the reddening listed in the E(B — V) column. 
Photometric errors are according to Table 1. Redshifts in boldface are spectroscopic. Candidates with no redshift are apparently hostless. The P(Ia) values of the candidates in the hostless and high-z samples 
were calculated assuming a uniform redshift probability distribution in the range < z < 0.35. /3 is the ecliptic latitude. 



fractions of SNe classified as Type la that we calculated for the 
four SN types. Each of the combinations is given a weight ac- 
cording to the number of permutations with the same fraction of 
SNe la. The outcome of this calculation is a probability distribu- 
tion for the real number of SNe in our sample, given that 16 SNe 
were classified as Type la. From this probability function we derive 
that the most probable "true" number of SNe la in our sample is 
Ma = 17+^ ±3.8, where the first error is the 68% ranse Poisson 
error, and the second one is a systematic 68% range error due to the 
uncertainty in classification. 



6.3 Supernova rate 

We now derive the luminosity-normalized SN la rate, ri a , in 
SNuband unit£|. The rate is calculated using 



la : 



Ma 



%Jm{z)j^{z)dv 



(2) 



where A^ a is the number of SNe la, dV is a comoving volume ele- 
ment, 7ii m is the luminosity density originating from galaxies which 
are above the SDSS limiting magnitude, and 77^ (z) is the effective 
visibility time (or "control time") of the z'-th image set, i.e., the time 
during which the SN is detectable. The integration is over the cos- 
mological volume in each set, and the summation is over image 



1 SNu w = SN (100 yr lO 10 !^)" 1 
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sets. Given e(m), the detection efficiency function as a function of 
magnitude m, 



(3) 



where, m e ff(f) is the effective SN light curve determined by the 
time difference, Afj, between the reference and registered images 
in each set i, 



meff (t) = -2.51og (lO" a4m W - W-°-Mt+^ . 



The mean redshift which we probe in this work, which de- 
pends on the visibility time and thus on the efficiency function, is 
given by 



(5) 



In order to calculate the luminosity density, j(z), we again 
use the Blanton et al. (2003) galaxy luminosity function. We con- 
vert their luminosity function, which is given for the SDSS bands 
shifted to z = 0.1, back to the rest-frame SDSS bands. We also 
account for luminosity evolution using their luminosity evolution 
parameter Q, thus obtaining the luminosity density as a function of 
redshift. Integrating over the luminosity function up to the limiting 
magnitude at each redshift provides jii m (z). 

The SN la rate must also be corrected for host extinction. 
Riello & Patat (2005) performed Monte Carlo simulations in which 
they modeled the dust distribution in host galaxies and accounted 
for various bulge-to-disc ratios and total optical depths. They found 
that the factor, /, by which SN la rates need to be corrected, is 
1 .27 < / < 1 .91, for Milky- Way-like dust. A similar, though lower, 
correction factor was derived by Neill et al. (2006), who derived 
the type la SN rate at z « 0.5. Neill et al. (2006) considered both 
Gaussian and exponential host extinction distributions in their de- 
tection efficiency simulations. They found a correction factor of 
1.10 < f < 1.37. Based on these studies, we adopt an intermediate 
correction factor of / = 1.25 to our SN rates. 

The derived SN la rates in the g and r bands, using Eq. 2 and 
also correcting for host extinction, are 

4 = (l 1.5+^+1-1 ± 2.5) x 10- 2 rf SNu r , 

r L = (l4.0^1;4±2.5)xl0- 2 ^ SNu g 

at a mean redshift of (z) = 0.20. The first error is due to the Poisson 
fluctuations in the SN number. The second is a systematic error due 
to the uncertainty in the detection efficiency function, calculated by 
using the efficiency function upper and lower (la) limits. The third 
error is the systematic classification error. 

For comparison with previously published rates, we convert 
our luminosity-normalized rates also to a volumetric rate. We do so 
by replacing Eq. 2 with 



Ma 



(6) 



./total (z) 



-dV 



where j tota \ is the total luminosity density. The resultant volumetric 
rates are then 



*L(0-2) 
*L(0-2) 



1.75 



+0.40+0.17 
-0.32-0.14 



1.89 



+0.42+0, 
-0.34-0, 



0.40^ X 10" 

}f±0.42) x 10" 



yr h^ Mpc 
yr _1 hj Mpc 



0.00015 



£0.0001 



a 



(4) S 



5xl0- 5 



oCappellaro et al. (1999) 
□ Madgwick et al. (2003) 
ABIanc et al. (2004) 
OHardin et al. (2000) 
xTonry et al. (2003) 
• Neill et al. (2006) 
<Barris & Tonry (2006) 
*Dahlen et al. (2004) 
>Botticella et al. (2008) 
oNeill et al. (2007) 
vDilday et al. (2007) 
■ This work 



■i 



5 



0.1 



0.2 



0.3 



0.4 



0.5 



Figure 7. Type la SN rates from different authors (see legend). Error bars 
are from quadrature additions of the various errors listed in Table 3. 



The difference between these two volumetric rates is due to the 
different evolution of the luminosity density in each band. This 
difference is an inherent weakness of deriving volumetric rates 
from luminosity-normalized rates, but in our case the difference 
is smaller than any of the other sources of uncertainty. 



7 COMPARISON WITH PREVIOUS MEASUREMENTS 

In this section, we compare our rate measurements to previously 
reported low-redshift SN rates. Most of these measurements (e.g., 
Cappellaro et al. 1999; Hardin et al. 2000; Blanc et al. 2004) were 
given in 5-band SNu units. These rates were then converted to 
volumetric rates using the luminosity density at the relevant red- 
shift. However, various luminosity functions were used to convert 
to volumetric rates. For example, Blanc et al. converted their rate 
and the rates of Cappellaro et al. (1999), Hardin et al. (2000), and 
Madgwick et al. (2003), using the 2dF redshift survey luminos- 
ity density (Cross et al. 2001). In contrast, Botticella et al. (2008) 
fit a set of luminosity density measurements (Norberg et al. 2002; 
Bell et al. 2004; Blanton et al. 2003; Faber et al. 2007; Tresse et 
al. 2007), with a smooth function of redshift, and used it to per- 
form the conversion to volumetric rates. We will repeat the Bot- 
ticella et al. (2008) conversion of volumetric rates of previously 
published luminosity-normalized rates (see Table 3), using their 
redshift-dependent luminosity density, 

j B (z) = (1.03 + 1.76 xz) x 10 8 Lq Mpc" 3 . (7) 

Dilday et al. (2008) have recently reported a Type la SN rate 
from the SDSS-II Supernova Survey (Frieman et al. 2008). The 
majority of SNe found in this survey, in contrast to our survey, 
have been confirmed spectroscopically (resulting in a lower red- 
shift range being probed). The SN la rate measured by Dilday et 
al., based on 17 SNe at z ~ 0.09 is higher than our rate measure- 
ment by a factor of ~ 1.5 but consistent within the errors. Table 
3 shows these various low-redshift rate measurements and Fig. 7 
shows a compilation of rate measurements to z < 0.5. 

Our rate measurement is consistent with other low-redshift 
rate measurements, and ranks with the most accurate among them. 
However, as mentioned above, some of the previously published 
volumetric rates, which appear in Table 3, would change, depend- 
ing on the luminosity density that is used to derive them. A remain- 
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Table 3. Comparison of low-redshift Type la SN rate measurements. 



< z > 


Ni a 




Rla 




Author 






h 2 10 SNu 5 


h 2 10 SNu, 


10~ 5 yr- 1 Mpc" 3 




~0 


70 


0.18±0.05 




1.85±0.5 


Cappellaro et al. (1999)* 


0.09 


17 




n 935+ - 07 
u.zjd_ 006 


9 q+0.9 
Z -^-0.7 


Dilday et al. (2008) a 


0.098 


19 


0.196±0.098 




2.4±1.2 


Madgwick et al. (2003)* 


0.13 


14 


n ioc+0.044+0.028 
u - iZJ -0.034-0.028 




i ^ro+0.56+0.35 
i - JO -0.43-0.35 


Blanc et al. (2004)* 


0.14 


4 


n 99+O.I7+O.O6 
u - zz -0. 10-0.03 




9 O+2.2+0.7 
z -°-l. 3-0.4 


Hardin et al. (2000)* 


0.2 


17 


14 +0.03+0.01 , 03 


14 +0.03+0.01 , 03 
UA ^-0.03-0.01 iUU ' 3 


1 fiO+0 .42+0.18 , ft 42 


This work 


0.25 


1 






1.7±1.7 


Barris & Tonry (2006) 



Notes: a Luminosity-normalized rate derived from a volumetric rate. * These rates have been converted to volumetric rates using the redshift-dependent 

luminosity density function from Botticella et al. (2008). 



ing uncertainty in such comparison arises from the fact that the 
rates in SNu of Cappellaro et al. (1999), Madgwick et al. (2003), 
Blanc et al. (2003), and Hardin et al. (2000) are multiplied by the 
luminosity density at the mean redshift of each survey, thus not tak- 
ing into account that the effective volume of the survey is a function 
of redshift. 

In order to compare our luminosity-normalized rate in SNug 
units also to rates given in SNug units, we adopt the Lupton et al. 
(2005) conversion between B and g band magnitudes, 

£ = g + 0.2271+0.313(g-r). (8) 

Using the mean luminosity- weighted color, g — r = 0.53, of the 
SDSS galaxy subset sample described in §4, together with Eq. 8, 
implies a factor of 1.03 increase going from SNug units to SNu# 
units. Using this factor, our rate measurement in SNug units is n a = 
(l4.4+|6+j-4 ±2 .6) x 10- 2 rf Q SNu fi . As seen in Table 3, this 
again agrees with previous measurements. 



tion efficiency function will remain the same. Although SNe found 
by archival search methods, such as ours, can not be studied spec- 
troscopically, a full SDSS -I sample could be useful for improved 
investigations of SN rates as a function of galaxy type and environ- 
ment. 
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8 SUMMARY 

We have conducted a low-redshift photometric SN survey using 
archival data from SDSS -I overlapping fields. Based on the number 
of Type la SNe that we find, Af Ia = 17+^ ±3.8, and keeping track 
of the various sources of error and bias, we have derived a SN la 
rate of rf a = (l4.0+ 2 J+}{ ±2.5) x 10~ 2 /z 2 SNu g , or a volumet- 
ric rate of (i-^o.'S-ais ± - 42 ) x 10 ~ 5 Y r_1 h io M Pc~ 3 - 
Our derived rates are consistent with previous measurements, but 
rank with the most accurate ones. However, SN la rates at low- 
redshift, including ours, still suffer from several sources of uncer- 
tainty. In our case, the uncertainty is mostly due to small numbers. 
The derivation of volumetric rates using different luminosity func- 
tions and different extinction corrections is another source of ambi- 
guity when comparing different measurements at similar redshifts, 
and when comparing observations with model predictions. Nev- 
ertheless, we have shown that there is a vast amount of archival 
SDSS data that can be used for studying SNe at a low cost. The 
full SDSS -I SN sample, once mined, would include several hun- 
dreds of SNe, comparable to the ~ 500 expected from SDSS-II, (a 
survey designed specifically for finding SNe). Assuming a similar 
fraction of Type la SNe as we found, both the Poisson and bino- 
mial classification uncertainties for such a large sample would be 
reduced to the ~ 5% level, while the uncertainty due to the detec- 
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